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ABSTRACT 

As part of the Transit Ephemeris Refinement and Monitoring Survey (TERMS), we present new 
radial velocities and photometry of the HD 192263 system. Our analysis of the already available Keck- 
HIRES and CORALIE radial velocity measurements together with the five new Keck measurements 
we report in this paper results in improved orbital parameters for the system. We derive constraints 
on the size and phase location of the transit window for HD 192263b, a Jupiter-mass planet with 
a period of 24.3587 ± 0.0022 days. We use 10 years of Automated Photoelectric Telescope (APT) 
photometry to analyze the stellar variability and search for planetary transits. We find continuing 
evidence of spot activity with periods near 23.4 days. The shape of the corresponding photometric 
variations changes over time, giving rise to not one but several Fourier peaks near this value. However, 
none of these frequencies coincides with the planet's orbital period and thus we find no evidence of 
star-planet interactions in the system. We attribute the ~23-day variability to stellar rotation. There 
are also indications of spot variations on longer (8 years) timescales. Finally, we use the photometric 
data to exclude transits for a planet with the predicted radius of 1.09 Rj, and as small as 0.79 Rj. 
Subject headings: planetary systems - techniques: photometric - techniques: radial velocities - stars: 
individual (HD 192263) - stars: activity - stars: spots 



1. INTRODUCTION 

Detecting exoplanets via the radial velocity (RV) 
method and subsequently monitoring their transit win- 
dows is the most fruitful strategy for exoplanets orbiting 
bright stars, which represent the best candidates for at- 
mospheric studies. Ground and space-based transit sur- 
veys have revealed nearly 200 transiting exoplanets, but 
most of those orbit faint stars, while the RV technique 
is best suited for brighter stars. Moreover, searching for 
transits of known RV planets allows the selective moni- 
toring of intermediate and long-period planets, of which 
only a few are known to transit so far. The Transit 
Ephemeris Refine ment and Monitoring Survey (TERMS; 
IKane et al.ll2009h aims to improve the orbital parameters 
of RV exoplanets and monitor them photometrically dur- 
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ing t heir thusly constrained transit windows. 

In Drag ornir et alj (|201lD . we presented new transits 
and improved parameters for several known transiting 
exop lanets, originally disc overed by the Super WASP sur- 
vey (jPollacco et al.l l2006h . We demonstrated that the 
photometric precision required to detect and characterize 
transits of giant planets is easily attainable by modest- 
sized, ground-based facilities such as the Cerro Tololo 
Inter- American Observatory (CTIO ) 1.0 m telescope. 
Through TERMS, the ephemerides of HD 156846b, HD 
114762b, HD 63454b and HD 168443b have been re- 
fined and transit searches conducted in each case (see 
IKane et al.ll2011al. IKane et alJl20lTH IKane et al.l 1201 Id 
and iPilvavskv et al.ll201ll respectively) . 

Exoplanets discovered using the RV method can some- 
times be controversial, and HD 1922 63b is one such mis- 
chief. It was first published in 2000 (jSantos et all 1200011 
as a planet with a period of 24.13 days and msini = 0.73 
Mj. These results arose from an analysis of RV measure- 
ments o btained using the CORALIE spectrograph. A 
paper by iVogt et all (|2000l ) followed, reporting a similar 
solution based on their Keck measurements, but noting 
that the chromospheric activity appeared to vary with 
a period cl ose to that of the suspected planet. Two 
years later, I Henry et al.l (|2002f ) attributed the RV sig- 
nal at least partly to stellar variability, as indicated by 
their photometric and spectrophotometric data. Indeed, 
a modulation with a period of ^24 days is clearly visible 
in the light curves, and the power spectrum of the Ca 
II H and K spectrophotometric observations exhibits a 
significant peak at the same period. In the end, un- 
like other plan et-like RV signals that were caused by 
stellar activity (lOueloz et al.ll200H; iDesidera et al.ll2004l: 
iPaulsonet atl I2004D. HP 192263b made a convincing 
comeback in iSantos et~aTI (|2003D . New CORALIE mea- 
surements demonstrated that the RV variation remained 
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coherent in amplitude and phase for over three years, 
while new photometric observations from La Palma re - 
vealed significant changes over time (jSantos et al.ll2003[ ). 

In this paper, we present new Keck RV observations 
which we use to refine the orbital parameters (section 
2) and the transit ephemeris of the revived HD 192263b 
(section 3). We introduce new photometry of the host 
star obtained between 2002 and 2011 in section 4, and 
report on the stellar variability during this period in sec- 
tion 5. Finally, in section 6 we analyze the photomet- 
ric measurements acquired during the predicted transit 
window and exclude transits with the predicted depth of 
2.5% with high confidence. We conclude in section 7. 

2. RV MEASUREMENTS AND REVISED ORBITAL MODEL 

2.1. Stellar Properties 

A K2 dwarf, HD 192263 is a relatively br ight star (V 
= 7.8 ) located at a distance of 19.3 pa rsecs (iKoen et al 
2010). It is also a BY Dra variable (|Kazarovets et al 



TABLE 1 
Stellar Properties 



2006), a class of active stars that experience brightness 



variations due to their spotted surface. 

We used Spectroscopy Made Easy (SME) on a spec- 
trum of HD 192263 (taken without th e iodine cell 
used f or RV measurements), as detailed in lValenti et al.l 
(120091). This pro c edure is based on that described in 
iValenti fc Fischer] ()2005[ ). with an added improvement 
for self-consistency between values of stellar properties 
obtained fr om spectroscopy an d those determined from 
isochrones (|Valenti et al.l 120091 ) . The improved method 
derives a value for the stellar surface gravity (log g), 
effective temperature (T e //), iron abundance ([Fe/H]) 
and alpha-element enrichment ([a/Fe]) from the spec- 
troscopic analysis. The last three parameters and the 
luminosity (L) are used to obtain a value for log g from 
isochrone models (log <foo)- The spectroscopic analysis 
is run again, now with log g fixed at the value of log gi SO 
from the previous iteration. The loop continues until log 
g and log gi SO agree. For HD 192263, the effective tem- 
perature resulting from this method is 4996 ± 44 K, and 
the stellar r adius (calculated fro m L and T e ft) is 0.73 
± 0.01 R©. iTakeda eTa l. (2007) carried out a Bayesian 
analysis using the stellar parameters arising directly from 
the SME procedure (log gi SO ,T e ff and [Fe/H]) and a 
dense grid of theoretical evolutionary tracks, and found 
a stellar radius of 0.77 ± 0.02 R . This value agrees 
with our SME result within t he 2a uncertainties. Fi- 
nally, Ivan Belle fc von Braunl ((2009) estimated the ra- 
dius of HD 192263 using spectral energy distribution 
(SED) fitting and obtained a value of 0.76 ± 0.02 R Q , 
which agrees with our SME result within the ler uncer- 
tainties. We adopt the SME stellar radius value for our 
calculations but compute a predicted transit dep th using 
the largest published value (jTakeda et al.ll200"7T ) as well, 
to show that transits can be excluded for both cases. 

From the SME analysis we also find values for the effec- 
tive temperature, surface gravity, iron abundance, pro- 
jected rotational velocity, and mass of the host star. The 
stellar properties are listed in Table 1. 

2.2. RV Measurements 

A total of 181 CORALIE (an echelle spectrograph 
on the l- 2 m Sw iss telescope at La Silla, Chile; 
lUdrv et al.l (I2000f0 RV me asurements have been re- 
ported in lSantos et al.l (|2003l ). of which those included in 



Parameter 



Value 



Reference 



V 
B-V 
Distance (pc) 
T efi (K) 
log 9 
[Fe/H] _ 
vs'mi (kms 1 ) 
logR' HK 
(M Q ) 
R* (Rq) 

age (Gyr) 



7.767 
0.957 
19.3 
4996 ± 44 
4.628 ±0.060 
0.054 ± 0.030 
< 1.0 ±0.50 

-4.56 
0.807 ± 0.015 
0.73 ±0.01 
2.09 ± 2.54 



Koen et al. (2010) 
Koen et al. (2010) 
Koen et al. (2010) 
This work 
This work 
This work 
This work 
Henry et al. (2002) 
This work 
This work 
This work 



TABLE 2 
CORALIE Radial Velocities 



Date 


Radial Velocity 


Uncertainty 


(BJD - 2440000) 


(ms -1 ) 


(m s — 1 ) 


11355.814795 


-64 


4 


11367.830021 


10 


12 


11375.760201 


-10 


13 


11381.750648 


-69 


7 


11384.720313 


-81 


6 


11390.738894 


16 


7 



Note. — Table 2 is presented in its entirety in the elec- 
tronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 



ISantos et all (pOOOl ) are a subset. In addition there are 
31 published Ke ck measurements ([Butler et al.l (|2006f ): 
iVogt et al.l (f2000h ). To these we add 5 new Keck observa- 
tions, acquired between 2006 and 2011. All Keck obser- 
vations were mad e with the HIRES echelle spectrograph 
(|Vogt et all I1994D on the 10-m Keck I telescope. The 
HIRES instrument uses an iodine cell through which the 
starlight passes before reaching the slit. Minute Doppler 
shifts in the features of the stellar spectrum are accu- 
rately measured against a wavelength reference provided 
by the dense set of absorption lines in the iodine spec- 
trum. The RVs were extracted f ollowing the procedure 
described in lHoward et al.1 (|2009l ). 

An offset of 10672 m s" 1 (the median of the CORALIE 
RVs) was added to each of the CORALIE RVs to place 
them on the same scale as the Keck RVs. The complete 
set of CORALIE and Keck RV measurements are given 
in Tables 2 and 3, respectively. We emphasize that all 
CORALIE measurements and all Keck measurements ex- 
cept the 5 we report in this paper, have been previously 
published. 

2.3. Keplerian Model 

We fit the set of 217 available RVs using a single-planet 
Keplerian model ba sed on the techniques described in 
IHoward et al.l ((20100 and the partially linearized, least- 
squares fitting method of I Wright fc Howa rd (2009;). An- 
alyzing the data from the two telescopes together rather 
than separately provides smaller uncertainties on the pe- 
riod and mid-transit time. We allowed for an RV offset 
between the Keck and CORALIE measurements (4.10 ± 
1.93 m s _1 ). We also allowed for an offset between Keck 
measurements taken before and after JD = 2453237, due 
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Keck Radial Velocities 




Date 


Radial Velocity 


Uncertainty 


(BJD - 2440000) 


(ms- 1 ) 


(ms -1 ) 


10984.063121 


32.75 


1.30 


11011.914062 


-11.45 


1.28 


11050.879437 


26.96 


1.38 


11069.890238 


-32.85 


1.16 


11312.086411 


-56.44 


1.23 


11313.112603 


-49.31 


1.30 


11342.057784 


16.00 


1.71 


11342.983844 


27.90 


2.17 


11367.915954 


39.38 


1.77 


11409.933630 


-39.58 


1.96 


11411.879482 


-13.05 


1.58 


11438.769456 


18.80 


1.60 


11439.828250 


28.04 


1.72 


11440.887793 


30.68 


1.87 


11441.833444 


42.11 


1.62 


11704.020541 


-30.66 


1.67 


11793.833579 


-55.58 


1.84 


12004.133959 


34.71 


2.13 


12008.138984 


7.54 


1.45 


12031.056206 


15.41 


1.54 


12063.060897 


-46.80 


1.72 


12094.907368 


-2.04 


1.54 


12128.910768 


30.51 


2.13 


12391.141713 


47.72 


1.87 


12536.755506 


31.62 


1.53 


12778.109319 


10.06 


1.36 


12833.916425 


42.13 


1.63 


12853.999751 


17.87 


1.59 


13181.024857 


-37.53 


1.33 


13239.861434 


-4.41 


1.04 


13546.990720 


-45.96 


1.05 


13969.022212 


-19.25 


0.83 


14810.722283 


1.86 


1.14 


15043.798793 


-5.86 


1.04 


15412.021408 


37.22 


0.85 


15782.862190 


33.76 


0.98 


Note. — The values in bold are the five 


new RV mea- 



surements we report in this paper. 



to a CCD and optics upgrade on that date (1.61 ± 5.84 
m s -1 ). The rms of the RV residuals is 13.14 m s -1 and 
the xled is 9 - 05 > likel y because HD 192263 is such an ac- 
tive star and its variability is not properly accounted for 
by the uncertainties on the RV measurements. To cor- 
rect for this, we added a jitter term (<7j) in quadrature to 
the measurement uncertainties (<Jrv)'- w 2 — <j 2 R V + a 2 , 
and u sed w instead of orv for the Keplerian fit ([Wright! 
2005). We chose a 3 = 10.15 m s to satisfy the con- 
dition x 2 — 1- Fur thermore, this value is equivalent to 
the excess found by lSantos et al.l (|2003| ) in the residuals 
of their RV measurements, and is within the range of 
10 - 30 m s _1 predicted by iSaar et al.l (|1998l ) from log 
R' HK = -4.56 (jHenrv et ai J 120021 ) . This added measure 
also reduces the rms to 12.42 m s _1 . 

The parameter uncertainties were determined from the 
sampling distribution of each parameter throug h a non- 
p arametric boo t strap analysis (|Freedmanlll981l) . 

ISantos et al.l (|2003| ) find a long term trend in their 
CORALIE data, to which they fit a line with a slope 
of 4.8 ± 0.8 m s" 1 yr" 1 . They are unable to determine 
the source of this trend, which could be due to the pres- 
ence of another companion, or a RV variation caused by 
stellar activity. We also include a linear velocity trend 
(dv/dt) in our models. The best-fit value for the offset 




-0.4 -0.2 0.0 0.2 0.4 

Orbital phase 

Fig. 1 — CORALIE (red diamonds) and Keck-HIRES (blue dots) 
radial velocities. The error bars shown are the original measure- 
ment uncertainties, with the jitter term (<trv) not included. For 
most Keck measurements, the error bars are smaller than the size of 
the data points. The best-fit orbital solution is overplotted (dashed 
line). The shaded region corresponds to the 3c transit window and 
phase 0.0 is the predicted time of mid-transit. See section 2 in the 
text for details. 

between the CORALIE and Keck RVs is now 5.14 ± 2.18 
m s _1 while the offset between Keck measurements be- 
fore and after JD = 2453237 becomes -17.67 ± 9.21 m 
s _1 . The best fit value for the linear velocity trend (when 
analyzing the Keck and CORALIE data combined) is 
0.0070±0.0017ms- 1 day -1 (or 2.56±0.62 m s^ 1 yr' 1 ). 
We note that including a trend does not significantly 
lower the \ 2 ■ When fitting only the CORALIE data we 
find a linear trend with a slope of 5. 21 ± 0.70 m s -1 yr" 1 
which agrees with that obtained by ISantos et al.l (|2003|) . 
However, fitting only the Keck data (which spans a longer 
time range encompassing that sampled by the CORALIE 
measurements) gives a slope of —0.37 ±0.73 m s _1 yr _1 , 
a value consistent with the absence of a trend. Finally, 
we repeated the analysis for each of the three data sets, 
this time using the original measurement uncertainties 
(no jitter term). The results are statistically consistent 
with those obtained with the jitter term included. We 
conclude that additional measurements are necessary to 
firmly determine the source of the linear trend present 
in the CORALIE data. For the reasons described above, 
and because the uncertainties on the resulting values for 
the period and mid-transit time are smaller, we adopt 
the solution without the trend. 

The parameters from our Keplerian fit (both with and 
without a linear trend as a free parameter) are given in 
Table 4, together with those reported by iButler et al.l 
(2006?) for comparison. The folded data and adopted 
model are plotted in Figure 1. Figure 2 shows a zoomed- 
in view of the shaded area (top panel) and the residuals 
to the RVs within this area (bottom panel). 

3. PREDICTED TRANSIT WINDOW AND 
CHARACTERISTICS 

From our newly derived stellar and planetary param- 
eters we can ascertain the properties of the pr e dicted 
transit. Using the models of iBodenheimer et al.l ([2003D 
and our revised measurement of the mass, we estimate 
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TABLE 4 
Keplerian Fit Parameters 



P ciXcin ic t c r 


This work - without trend 
(adopted) 


This work - with trend 


Butler et al. 2006 


P (days) 


24.3587 ± 0.0022 


24.3581 ± 0.0028 


24.3556 ± 0.0046 


T c a (BJD - 2440000) 


10986.74 ±0.21 


10986.73 ±0.21 


10987.22 ± 0.39 


T p b (BJD - 2440000) 


11796.9 ±6.8 


11795.0 ±4.6 


11994.3 ±3.9 


e 


0.008 ±0.014 


0.008 ± 0.014 


0.055 ± 0.039 


K (ms" 1 ) 


59.3 ± 1.2 


58.7 ±1.2 


51.9 ±2.6 


w (deg) 


184 ± 90 


157 ± 88 


200 


dv/dt (ms -1 jr" 1 ) 




2.56 ±0.62 




M p sin i (M;) 


0.733 ± 0.015 


0.726 ±0.014 


0.641 


a (AU) 


0.15312 ± 0.00095 


0.15312 ±0.00095 


0.15 


rms (ms ) 


12.42 


11.87 


12.5 



a Time of transit. 

b Time of periastron passage. 
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Fig. 2.— Top panel: CORALIE (red diamonds) and Keck-HIRES 
(blue dots) radial velocities between the edges of the 3cr transit win- 
dow (shaded region in Figure 1), with the best-fit orbital solution 
overplotted (dashed line). Bottom panel: Residuals to the radial 
velocities, over the same phase range (rms = 12.04 m s~ x ). The 
error bars shown are the original measurement uncertainties, with 
the jitter term (<trv) not included. 



a radius for the planet of R p = 1.09 R,/. It was shown 
bv lKane fc von Braunl (|20085 that the transit probability 
is a strong function of eccentricity and periastron argu- 
ment. The eccentricity of this orbit is small enough such 



that this introduces a minor effect. Based upon our pa- 
rameters (listed in the first column of Table 4) , the tran- 
sit probability is 2.49%, the predicted duration is 0.192 
days, and a predicted depth is 2.53%, easily within the 
range of our photometric precision (rms = 0.0065 within 
the 3<T transit window; see section 6 for details). 

We derived predicted transit times using the non- 
parametric bootstrap analysis described in Section 2.3. 
The predicted transit mid-point time used to fold the 
data was 2455882.84 ± 0.37 (JD). Note that the final 
Keck measurement obtained reduced the uncertainty on 
the transit mid-point from 0.47 days to 0.37 days, em- 
phasizing the importance of extending the time baseline 
of the RV measurements when attempting to reduce the 
total size of the transit window. In this case, the la tran- 
sit window has a total duration of 0.932 days. Although 
this is dominated by the transit mid-point uncertainty, it 
is small enough that we are able to attempt the detection 
of the transit. 

As mentioned in Section 2.1, iTakeda et~aTl (|2007t ) de- 
rive a slightly larger stellar radius of 0.77 ± 0.02 R Q . Us- 
ing this radius results in a transit probability of 2.61% 
and a transit depth of 2.28%. This is still well within 
our photometric precision such that we are able to make 
a definitive statement regarding the transit exclusion of 
this planet. 

4. PHOTOMETRY 
4.1. APT photometry 

We acquired 985 photometric observations using the 
Til 0.8 m automatic photometric telescope (APT), lo- 
cated at Fairborn Observatory in southern Arizona. The 
data span just over a decade, from April 13, 2001 to 
November 23, 2011. Measurements were obtained simul- 
taneously in the Stromgren b and y passbands by two 
EMI 9124QB photomultiplier tubes. The individual b 
and y differential magnitudes arc averaged to obtain the 
quantity A (6 + y)/2. The observing and data reduc- 
tion p rocedures are identical to those described in lHenrvl 
(fl999l) for the T8 0.8 m APT. 

The differential magnitudes A(b + y)/2 were converted 
to relative fluxes and normalized to 1 for the plots in 
Figure 3. Two comparison stars were considered for the 
differential photometry. Comparison star 1 (CI) is HD 
193328 (V = 7.48, B-V = 0.12, A2), while comparison 
star 2 (C2) is HD 193225 (V = 7.35, B-V = 0.29, F0). 
Typical precision of a single relative flux measurement 
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Fig. 3.— Til APT observations of HD 192263 covering 11 ob- 
serving seasons from April 2001 to November 2011. The top panel: 
the photometry for HD 192263, obtained using comparison star 
CI (rms = 0.0085). Middle panel: same as in the top panel, 
but with each half observing season normalized to the same mean 
(rms = 0.0067). Bottom panel: photometry of comparison star 
CI relative to comparison star C2 (rms = 0.0023). It is clear that 
the large scatter in the HD 192263 observations is mainly due to 
the target star itself, not to any variability in the comparison stars. 
The la error bar is shown in each panel. See sections 4 and 5 in 
the text for details. 



from Til is 0.0010 - 0.0020, depending on the quality 
of the night and the air mass of the observations. The 
standard deviation of the C1/C2 relative fluxes (see Fig- 
ure 3, bottom panel) is slightly larger than this (0.0023). 
This is probably a combination of the air mass, since 
HD 192263 lies near the celestial equator, and also sus- 
pected low-amplitude variability in C2. For this reason, 
we analyzed the differential photometry of HD 192263 
with respect to CI (see Figure 3, top panel). The 985 
A(6 + y)/2 magnitudes for HD 192263-C1 and C1-C2 are 
listed in Table 5. 
HD 192263 is observable from Arizona only during part 



TABLE 5 

Photometric Observations of HD 192263 from the Til 
APT 



Heliocentric Julian Date 
(HJD - 2,400,000) 


(HD 192263-Cl)i,j, 
(mag) 


(CI - C2) by 
(mag) 


52,012.9927 


0.5293 


0.0729 


52,020.9783 


0.5337 


0.0765 


52,022.9700 


0.5368 


0.0710 


52,023.9680 


0.5260 


0.0735 


52,025.9669 


0.5191 


0.0750 


52,027.9535 


0.5215 


0.0723 



Note. — Table 5 is presented in its entirety in the elec- 
tronic edition of the Astrophysical Journal. A portion is 
shown here for guidance regarding its form and content. 

of the year, which explains the larger gaps in the data 
between observing seasons. In addition, the star is at 
opposition around July 24, during the annual Summer 
Shutdown of the APTs. This gives rise to the shorter (8 
- 10 weeks) gaps between the first and second clusters of 
data points (half seasons) , between the third and fourth, 
etc. As such, the first two sections of the light curve 
correspond to the first observing season, the third and 
fourth belong to the second observing season, and so on. 

We prepared the data for the transit search by normal- 
izing it such that each half season has the same mean 
value (see Figure 3, middle panel). This normalization 
affects the data on a timescale of months, thus preserv- 
ing the shape and depth of any potential transits which 
would last less than 8 hours (see section 3). 

4.2. ASAS photometry 

HD 192263 was also observ ed as part of the All Sky 
Automated Surve £3 (ASAS) (|Poimanskilll997l) . A total 
of 345 data points are available in the ASAS-3 photo- 
metric V-band catalogue. These observations were col- 
lected between March 27, 2001 and November 10, 2009. 
The quality and cadence of this data set are lower than 
for the APT photometry. As a consequence, we use the 
ASAS observations to test our conclusions regarding the 
variability of the star, but we do not employ them for 
the transit search. 

5. VARIABILITY OF THE HOST STAR 

In this section, we describe an investigation of the stel- 
lar variability of HD 192263 based on the 11 seasons 
of APT photometry described in the previous section. 
We generate an amplitude spectrunri from a discrete 
Fourier transform of the time series (see top panel of 
Figure 4) and search it for significant peaks. Beyond 
the frequency range shown in Figure 4 (0.00 to 0.07 cy- 
cles day -1 ), only harmonics of the significant frequencies 
shown in the top plot rise above the noise. 

The middle panel and the inset of Figure 4 show the 
spectral window functions for the two significant frequen- 
cies shown in the top panel. These indicate the aliases 
introduced into the amplitude spectrum by the observ- 
ing cadence, with regular gaps on nightly, biannual and 
annual timescales. 

11 http:/ /www.astrouw.edu.pl/asas/?page=aasc&catsrc=asas3 

12 An amplitude spectrum shows the amplitude of signals present 
in the data plotted versus frequency. A power spectrum is obtained 
by plotting the square of the amplitude versus frequency. 
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Fig. 4. — Fourier amplitude spectrum (top), window functions for 
the two significant frequencies described in the text (middle and 
inset) and Fourier spectrum of the light curve residuals after pre- 
whitening (bottom) of the APT photometry. The dotted, dashed 
and solid horizontal lines in the top and bottom panels correspond 
to false-alarm probabilities of 4.6%, .3% and 0.1% , respectively, 
calculated using the SigSpec routine (Rccgcn 2007). Details can 
be found in section 5. 

The strongest peak in the top panel of Figure 4 has a 
frequency of 0.0003354 ± 0.0000055 cycles day" 1 , which 
is equivalent to a period of 8.17 years (or 2982 days ± 
50). This signal is clearly visible in the photometry, as 
shown in the top panel of Figure 3. In light of the already 



known BY Dra nature of the star, this 8.2-year timescale 
is consistent with an activity cycle in an early K dwarf 
like HD 192263, but its cyclical or periodic nature can 
only be determined with additional observations. The 
brightness changes are the result of changes in the filling 
factor of spots or active regions on the stellar photosphere 
(lOzdarcan et al.ll2010t)- Th ey may be driven by a stellar 



dynamo (Strassmcicr 2005), or may arise from random 
variations in a large number of spots on the surface of an 
active star without the requirement of a driving mecha- 
nism, as shown by the simulations of Eat on et al.l {1996). 

The amplitude spectrum further shows a cluster of 
peaks between 0.03 and 0.06 cycles day -1 , of which the 
largest occurs at 0.0427475 ± 0.0000084 cycles day" 1 
(period = 23.3932 ± 0.0046 days) with an amplitude of 
0.005 mag. The window function for a sinusoid of this 
frequency and amplitude (see inset of Figure 4) reveals 
that not all of these peaks are aliases, but that there 
are other independent frequencies present in this range. 
We first " pre- whitened" the light curve by sequentially 
removing the frequencies of the largest peaks (and their 
harmonics), until we reached the noise level. Six frequen- 
cies in this range were identified; they are listed in Ta- 
ble 6, along with the very low frequency discussed above. 
We then performed simultaneous least-squares fits to the 
data, where the frequencies, amplitudes and phases of 
these peaks were allowed to float, with the values ob- 
tained from the amplitude spectrum as initial guesses. 
The Fourier amplitude spectrum of the light curve resid- 
uals following the pre-whitening is plotted in the bottom 
panel of Figure 4. None of the identified frequencies cor- 
responds to the orbital frequency of the exoplanet HD 
192263 b. In fact, when we include that frequency and 
allow the least-squares fit to iterate, the frequency mi- 
grates to one of those listed in Table 6. Despite the fact 
that there are at least six frequencies in a narrow range, 
it is clear that the orbital frequency of the planet is not 
one of them. 

The Fourier peaks clustered around a period of about 
23.4 days are a result of the evolving nature of star spots. 
As spots or groups of spots form, evolve and eventually 
disappear, contemporaneously or successively, the shape 
and amplitude of the light curve change over the course 
of one to a few stellar rotation periods. Thus, although 
the true stellar rotation period is close to 23.4 days, it 
cannot be determined to the frequency resolution of the 
time series. 

The frequency analysis and the behaviour of the light 
curve are consistent with rotational modulation of the 
light output of an active KO dwarf observed for about 
a decade. If the star is undergoing an activity cycle, 
then the spottedness will also change on that timescale 
(varying the amplitude of rotational modulation). In the 
light curve, shown in the upper two panels of Figure 3, 
one can see the changing width of the envelope of points 
over the 11-year time span. If the variation at a pe- 
riod around 23.4 days were constant in amplitude, that 
envelope would maintain the same width (since the pho- 
tometric scatter of the comparison star remains uniform 
throughout the data set, as can be seen in the bottom 
panel of Figure 3). All indications are that the ampli- 
tude of this signal changes over the span of the ^8-year 
variation, and not in a random way. 
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TABLE 6 

List of observed frequencies for. HD 192263 



# Frequency (c/d) Period (days) Amplitude (rel. flux) 





<j frequency 


(j period 


(7 cLITlplitU.dc 


1 


0.0003354 


2982 


0.00612 




0.0000055 


50 


0.00059 


2 


0.037823 


26.439 


0.00196 




0.000024 


0.017 


0.00034 


3 


0.040741 


24.5640 


0.00273 




0.000017 


0.0085 


0.00042 


1 


0.0424 


23.5849 


0.0024 




0.0015 


0.8650 


0.0012 


5 


0.0427475 


23.3932 


0.00278 




0.0000084 


0.0046 


0.00012 


6 


0.043625 


22.9226 


0.00333 




0.000017 


0.0090 


0.00040 


7 


0.04485 


22.30 


0.00256 




0.00029 


0.15 


0.00095 



Note. - The full least-squares fit solution including 
phases for these frequencies is available upon request. 



1.02 



1.01 
I 1 - 00 

(D 

0.99 



0.98 



-0.4 -0.2 0.0 0.2 0.4 

Rotational phase 




Our findings indicate that the periodic variability re- 
ported by Henry et al. (2002) and Santos et al. (2003) 
persists. To further verify this, we subtracted the 8.17- 
year signal from the data and phased the residuals at a 
period of 23.39 days (the largest amplitude signal in this 
period range), as shown in the top panel of Figure 5. 
A dominant variation, which remains roughly in phase 
throughout the time span of the photometry, stands out 
in this phase diagram. 

We also see periodicity around 12 days in some of the 
half-season data sets, and we see the first harmonics of 
the periods near 23.4 days in the Fourier spectrum. This 
is in agreement with the period of 12.2 ± 0.1 days ob- 
tained by Henry et al. (2002) from their spectrophoto- 
metric Ca II H and K observations. They proposed that 
the 12.2-day period is half the stellar rotation period, 
arising in the data when active regions are present on 
opposite hemispheres of the star. 

If we take a value of 23.4 days as the rotation period 
of the star and our stellar radius of 0.73 Rq, we find a 
rotational velocity of 1.58 km/s, which is consistent with 
our measured value of vshii < 1.0 ± 0.5 km s _1 (i.e. 
v ro t > v sini). 

Using photometry spanning just over a year, Henry et 
al. (2002) found a rotation period of 24.5 ± 0.5 days, 
which agrees with the period of the planet within the 
uncertainties. This led Santos et al. (2003) to suggest 
the possibility of star-exoplanet interactions in the sys- 
tem. However, none of our new, more precise values of 
the dominant frequencies associated with rotation, based 
on a much longer data set, match the orbital period of 
the planet. Our least-squares test, trying to force a fit 
including the orbital period of the planet (24.3587 days), 
show that it is not part of a valid solution for the data. 
As an additional reality check, we phased the photom- 
etry at this period and found that the coherent signal 
visible in the top panel of Figure 5 disappears, as can 
be seen in the bottom panel of the same figure. Based 
on our data, we conclude that there is no indication of 
star-planet interactions in the HD 192263 system. 

For completeness, we have also examined the ASAS 
photometry with the goal of verifying the conclusions 
described above. This time series has a lower cadence, 
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Fig. 5. — APT photometry with 8.17-year trend removed, phased 
at the rotation period of the star (23.39 days; top panel) and at 
the orbital period of the planet (24.36 days; bottom panel). The 
coherent signal visible in the top panel disappears in the bottom 
panel, indicating that it is more likely to correspond to the stellar 
rotation period than to star-planet interaction. The la error bar 
is shown in each panel. Details can be found in section 5. 



covers a shorter time frame and its rms is three times 
larger than for the APT data. We find that while the 
two data sets are consistent with each other, the lower 
quality of the ASAS photometry prevents it from setting 
additional constraints on the variability of the star. 

6. TRANSIT EXCLUSION AND IMPLICATIONS 

The long APT photometric dataset for HD 192263 
makes possible the full coverage of the 3tr transit window. 
Figure 6 shows the photometry phased at the planet's or- 
bital period. Phase 0.0 is the location of the predicted 
mid-transit time of the planet. The solid line represents 
the predi cted transit signature, based on the analytic 
models of iMandel fc Agoll (j2002) and the values for the 
transit depth and duration calculated in section 3. The 
short dashed line corresponds to the predicted transit 
depth if the stellar radius is 0.77 R© (see section 3). At 
the top the entire orbital phase is shown, while the hori- 
zontal range of the bottom plot corresponds to the size of 
the 3<7 window. In the bottom panel, the vertical dashed 
and dotted lines indicate the extent of the la and 2a 
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transit window, respectively. 

The rms of the photometry in Figure 6 is 0.0065. We 
can therefore exclude edge-on (i = 90°) transits with the 
predicted depth of 0.0253 at the 3.9a level. (If we assume 
a stellar radius of 0.77 R Q , then the predicted depth 
is 0.0228 and we can exclude such transits at the 3.5a 
level.) While a grazing transit for this planet requires 
an orbital inclination of 88.7°, the cadence of the obser- 
vations only allows us to exclude transits corresponding 
to i >88.9°. This is equivalent to an impact parameter 
(6) < 0.86, and a transit duration > 0.01 days = 0.004 
orbital phase. The photometric precision is sufficient to 
rule out transiting planets with radii as small as 0.79 Rj 
at 2a confidence. The density of a planet with such a ra- 
di us would fall outside the ra nge predicted by the models 
of lBodenheimer et all (120031 ). 

The case of HD 192263b shows that deep and more 
importantly, long transits can be detected or ruled out 
using existing low-cadence photometry if the dataset cov- 
ers a sufficiently long time period. An additional benefit 
of a long time series is the full coverage of the 3a tran- 
sit window, providing a higher confidence in the transit 
exclusion in the case of a non-detection. The utility of 
the ephemeris refinement component of the TERMS ap- 
proach is clear: smaller uncertainties on the orbital pe- 
riod and transit time lead to a shorter transit window 
within which we need to assess whether the phase-folded 
photometry has adequate cadence and precision. 

We perform an additional test to check the presence or 
absence of a transit by comparing the predicted ampli- 
tude of the Rossiter-McLaughlin (RM) effect for this sys- 
tem with our RV measurements. For a system in which 
the orbital plane of the planet is al igned with the equator 
of the star, iGaudi fc Winnl (|2007| ) show that the ampli- 
tude of the RM effect is given by 



K R = vsmi- - (1) 

1 — 7 2 

where 7 = R p /R-,,. For 7 <C 1, equation 1 becomes 

K R = vwxi(^\ (2) 

Using a transit depth (R p /R+) 2 of 0.0253 and our up- 
per limit on vsini, we find Kr < 25.3 m s _1 . Given the 
rms of our RV residuals within the 3a transit window, 
we can exclude an RM effect with amplitude > 12.0 m 
s _1 . While these limits place a much looser constraint 
on the existence of a transit than the photometry does, 
they are consistent with the absence of a transit. 

7. CONCLUSIONS 

The planet orbiting HD 192263 is an example of the 
challenges posed by the detection of companions around 
active stars. In this paper, we presented five new RV 
measurements which we used to refine the orbital pa- 
rameters of HD 192263b. The new measurements pro- 
vide continuing support for the exi stence of the planet , 
in agreement with the conclusion of iSantos et al.l ([2003). 

We have also shown new photometry of the system. 
We perform a Fourier analysis of the photometry and 
find evidence for variability with periods near 23 days, 



1.02 



1.01 



ID 

0.99 



0.98 



-0.4 


-0.2 


0.0 
Orbital p 


iase 


0.2 


0.4 


* 

Z * 














: ¥ — * 








"* 


" * ■ • 




T , I 


* * 

liii 


I 


• • 


I \ 

i , 7 


-0.04 


-0.02 


0.00 


0.02 


0.04 



Orbital phase 

Fig. 6. — APT photometry, with each half-season normalized to 
the same mean (see text for details) and phased at the orbital pe- 
riod of the planet. The solid line is the predicted transit signature. 
The overplotted short dashed line is the predicted transit signa- 
ture if the star has a radius of 0.77 Rq. Top panel: full orbital 
phase. Bottom panel: the 3<r transit window (horizontal extent of 
the plot), with the vertical dashed and dotted lines enclosing the 
lcr and 2 a transit windows, respectively. The lcr error bar is shown 
in each panel. See section 6 in the text for details. 

which agrees with previous reports of stellar variability 
and which we attribute to stellar rotation. A detailed 
examination of the Fourier spectrum near this value re- 
veals a multiplet of peaks with periods ranging between 
22 and 27 days. The dominant signal in this cluster has 
a period of 23.3932 ± 0.0046 days. The identified sur- 
rounding peaks arise from the evolving nature of spots on 
the stellar surface, which affect the shape and amplitude 
of the light curve on the timescale of the stellar rotation 
period. Nevertheless, neither the dominant period nor 
any of the five other peaks match the orbital period of 
the planet (P = 24.3587 ± 0.0022 days), so we find no 
evidence of star-planet interactions. 

We also observe a longer-term trend which may be a 
~ 8- year activity cycle (if the cycle repeats), or may just 
be part of a longer interval of random variations in the 
star's spottedness. Continuing long-term monitoring of 
the star should more convincingly discriminate between 
the two scenarios. It is noteworthy that we do not see 
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a long-term variation in the RV measurements in phase 
with this long-term photometric trend. 

As our photometric dataset spans approximately a 
decade, we have good coverage of the 3<t transit window 
when the data are phased to the orbital period of the 
planet. Thus we are able to thoroughly exclude transits 
of the predicted depth (2.53%, corresponding to a planet 
with a radius of 1.09 Rj) for a planet with a mass of 
0.733 Mj. The absence of a detectable (> 12.0 m s" 1 ) 
Rossiter-McLaughlin effect is consistent with this result. 
We also exclude transit depths as low as 1.3% (corre- 
sponding to a planetary radius of 0.79 Rj). In the case 
of a non-edge-on orbital configuration, the cadence of the 
data allows us to rule out transits with impact parameter 
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